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second order statistics for 
quadrature, branches of a 

received transmitted siganJ as input to the receiver, a coarse synchronization and a coarse frequency offset compensation have being 
performed on the signal. It comprises a means for derotation, means for separation, means for filtering, means for estimating and 
means for detecting transmitted symbols in the received signal. The invention thereby improving co-channel rejection in wireless 
communication, thus making it possible to increase the number of communication channels for frequencies used. 



(57) Abstract: The present invention relates to a method for a digital receiver and a receiver exploiting j 
adaptive co-channel interference rejection in wireless communication. It uses digitally I, in phase, and Q, 
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A R eceiv er 
TECHNICAL FIELD 

Thepresent invention relates to a method and a receiver for exploiting second order 

statistics for adaptive co-channel interference rejection in wireless communication. 

5 DESCRIPTION OF THE BACKGROUND ART 

Co-channel interference limits the number of simultaneous users in today's wireless 
systems; techniques for efficient co-channel interference rejection in wireless communication 
systems are of great importance. Recently, antenna arrays have been proposed to obtain extra . 
degrees of freedom, so that spatio-temporal processing may be used for interference rejection. In 

10 some cases it is infeasible or unattractive to employ antenna arrays. 

The increased number of users in wireless systems motivates the research in the area of 
interference rejection techniques. In order to suppress interfering co-channel signals some degree 
of freedom in the received signal is generally exploited This freedom may be obtained by using 
multiple receiving antennas at the base station and/or by exploiting temporal structure in the 

1 5 communication signal. 

Use of antenna arrays in wireless systems has been intensively investigated recently; 
see P. Zetterberg and B. Ottersten, 'The spectrum efficiency of a basestation antenna array 
system for spatially selective transmission", IEEE Transactions on Vehicular Technology, 44(3): 
65 1-660, August 1 995, F.Pipon, P. Chevalier, P. Vila and J. -J. Monot, "Joint spatial and 

20 temporal equalisation for channels with ISI - theoretical and experimental results for a base 

i 

station antenna array", in Conference proceedings of Signal Processing Advances in Wireless 
Communications (SPAWC), pages 309-3 12, Paris, April 1997 IEEE, or D. Asztely and B. 
Ottersten, " MLSE and spatio-temporal interference rejection combining with antenna arrays", in 
DC European Signal processing Conference (EUSIPCO-98), island of Rhodes, Greece, 

25 September 8-11 1998. 

Multi-dimensional antenna array signal processing is traditionally focused on utilising 
the spatial dimension present in array output data. The temporal dimension (delay-spread) has 
until recently been dealt with in separate circuits (channel equalisation). Combined temporal and 
spatial processing, as noted in D. T. M. Slock, "An interference cancelling multichannel matched 

30 filter", in Communications Theosy Mini-Conference in conjunction with Globecom, pages 214 - 
218, New Yoric, NY, 1996, IEEE and in A J. Paulray and C B. Papadias, "Space-time 
processing for wireless communications", IEEE Signal Processing Magazine, 14(6): 49-83, 
1997, is advantageous when constructing highly efficient receiver methods for interference 
rejection. Drawbacks with antenna arrays are the increase in hardware and the difficulty to 
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implement them on mobile stations. The present invention describes a method and a receiver for 
adaptive co-channel interference rejection using one or more antennas. 

The method and the receiver according to the present invention are set forth through the 
attached independent claims. Further embodiments of the invention are set forth through attached 
5 dependent claims. 



SUMMARY OF THE INVENTION 

The present invention relates to a method and a receiver with a co-channel interference 
rejection scheme that exploits second order statistics for adaptive co-channel interference 

1 0 rejection in wireless communication. The invention claimed exploits the properties of certain 
signal constellations to obtain extra degrees of freedom that enable the application of powerful 
algorithms for interference rejection. This method may be applied to one or multiple coherent 
radio channels or one or more antennas. Examples of signal constellations include PAM, BPSK, 
MSK and also GMSKL Digitally I (in phase) and Q (quadrature) branches of a received 

1 5 transmitted signal are used as input to the receiver, a coarse synchronisation and a coarse 

frequency offset compensation have being performed on said signal A receiver according to the 
invention claimed comprises means for denotation, means for separation, means for filtering, 
means for estimating and means for detecting transmitted symbols in the received signal. 
Derotating the signal in said means for denotation outputs a complex valued time 

20 discrete base band representation of the received signal if a used modulation fonn is one 

dimensional. If a used modulation fonn is multiple dimensional, from possibly two vector valued 
signal sequences, a complex valued signal is separated in said means for separating it into its real 
and imaginary parts which vector is output from said means for separation containing the 
transmitted signal, co-channel interfering user signals, additive noise and other possible 

25 disturbance. 

The method according to the invention claimed comprises the following steps: 
estimating in said means for estimation, using the outputs from said means for separation, means 
for filtering, from said means for detecting as input, and performing one or more of the following 
quantities: fine synchronisation, fine frequency offset estimation compensation, or data model 
3 0 estimation, outputting one or more of said quantities to said means for derotating, filtering and 
detecting. 

Filtering in said means for filtering, uses the'real valued output from said means 

* • 

for separation and coefficients in said data model output from the means fof estimation and 
performing a whitening operation on the signal additionally taking into account the second 
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order properties of noise and co-channel interference, the filtered signal being the output to 
the means for estimation and detection. Said filtering means perform a whitening operation 
that takes into account the second order statistical properties of the noise and co-channel 
interference. 

5 Detecting transmitted symbols in said means for detecting, making use of the 

output from said means for filtering and the data model from the means for estimation, 
thereby improving co-channel rejection in wireless communication, thus making it possible to 
increase the number of frequencies used. 

In one embodiment of the present invention the model for the communication 

* 

10 channeHs a linear FIR channel. 

In another embodiment of the invention a model for interference and noise is a 
linear vector valued moving average process. In another embodiment a model for interference 
and noise is a vector valued autoregressive process. 

A further embodiment of the present invention sequence estimation is performed 
IS on the whitened sequence using a filtered channel and a residual covariance matrix. In another 
embodiment detection is performed by taking a binary decision, symbol by symbol based on a 
scalar input signal. In yet another embodiment detection is performed by a reduced 
complexity sequence estimation performed on the whitened sequence. In yet another 
embodiment linear equalization is performed followed by symbol by symbol detection, 

20 Data model parameters are in one embodiment of the invention estimated using 

training data. In another embodiment data model parameters are estimated using estimated 
symbols. 

■ 

In a still further embodiment of the invention the model order is prespecifed. In 
another embodiment the model order is chosen according to a model order selection criteria. 
25 In yet another embodiment the model order is adapted during communication according to 
error rate on the channel. * 

In yet another embodiment multi models are estimated and used in the filtering 

- 

and detection means, the model resulting in the lowest error rate providing the symbol 
estimates. 

■ * 

■ 

30 In a yet further embodiment of the invention timing estimation is performed on 











HI! 



training data. In another embodiment timing estimation may be performed using estimated 
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symbols. In yet another embodiment timing estimation is performed simultaneously as the 
data model parameter are estimated using training data. 

A further embodiment provides that frequency offset estimation is performed on 
training data. In another embodiment frequency offset estimation may be performed using 
5 estimated symbols. 

In yet another embodiment frequency offset estimation is performed simultaneously 
as the data model parameter are estimated using training data. 

In another embodiment frequency offset estimation is performed on training data. 
Frequency offset estimation may be performed using estimated symbols. 
10 In another embodiment frequency offset estimation is performed simultaneously as the 

datamodel parameters are estimated using training data. 

hi a still further embodiment said filtering means comprises a vector-valued linear IIR 
Infinite Impulse Response filter. Said filtering means performs a whitening operation that takes 

4 

into account the second order statistical properties of the noise and co-channel interference. 
1 5 Multichannel maximum-likelihood sequence estimation may be performed by applying a Viterbi 
algorithm to the whitened sequence. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 schematically illustrates a prior art receiver, and 
20 Fig. 2 schematically illustrates a receiver according to the present invention. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS OF THE INVENTION 
The present invention relates to a method and a receiver with co-channel 

■ « » * 

interference rejection scheme that exploits second order statistics for adaptive co-channel 
25 interference rejection in wireless communication. To illustrate the invention, a radio receiver for 
GSM is described. Extensions to other digital modulation schemes using signal constellations 
with the specified structure are included in the present invention. Communication signals are in 
general very structured, for example, the special temporal shape of a transmitted signal in Code 
Division Multiple Access schemes (CDMA), controlled by different user specific codes, used to 
30 separate concurrent channel users. 

The communication signal is in this case thus extremely useful to suppress 
unwanted signals. Depending on timing offsets, signal strength, etc, the cod#g separates users 
more or less efficiently. Herein, time division multiple access systems, are used to illustrated the 
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invention. For CDMA systems that use signal constellations with the specified structure (for 
example BPSK), receivers with adaptive interference rejection according to this invention can 
also be implemented. It is also shown how the structure of Binary Phase Shift Keying (BPSK), 
Pulse Amplitude Modulation (PAM), Minimum Shift Keying (MSK) and Gaussian Minimum 
Shift Keying (GMSK) communication signals can be exploited to increase the number of 
processing degrees of freedom. It is shown that, starting with a one d imensional symbol sampled 
lunication system, a two dimensional system may be obtained, see M. Kristensson, D. T. 



Slock and B. Ottersten, "Blind subspace identification of a BPSK communication channel", in 
Proceedings of the 30* Asilomar Conference on Signal, Systems and Computers, Pacific Grove, 

• * 

10 CA, November 1996, This facilitates the use of powerful spatio-temporal algorithms when 

detecting the user of interest; this even when there is only one available antenna and no temporal 

over-sampling is used. 

Fig. 1 illustrates a procedure that is well known and used in most modem digital 

receivers. Actual implementation and order can take many forms, the filtering and down 
15 conversion being accomplished in stages, sometimes digitally. The figure displays an overview 

of the first analogue blocks in a radio receiver. These blocks filter 14, 10, down convert and 

sample 12 a received signal r(/). 

In the present invention it is assumed that the digitally sampled J (in phase) and Q 

(quadrature) branches of the received signal, represented by ri(k) 16 and tq(^ 1 8, are available 
20 in the receiver. 

It is assumed according to the present invention that a coarse synchronisation and 
coarse frequency oflset compensation have been performed. An index, k, represents time and is 
assumed to be^^d^^^The quantities r f (jt) and r Q (Jfc) may be vector valued in case of 
multiple antennas and/or oversampling with respect to symbol time. The dimension ofthe 
25 vectors is equal to the product of the number of antennas and the oversampling factor. This is 
often referred to as the time discrete, digital base band communication signal The sampling rate 
in this first block in the receiver can be faster than the symbol rate. Oversampling with respect to 

» 

the symbol period is allowed 

Below, blocks in the digital receiver involved in the methods according to the 
30 present invention to exploit the modulation form for interference rejection are described. The 
following blocks are central to the invention and are depicted in Fig. 2. 
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Depending on the modulation for the transmitted signal, the Derotation block 20 
functions in one of two ways. 

First, if the transmitted modulation form is one dimensional, e.g. only needs one 

* * 

basis function to be described, then the output of the Derotation block 20 in Fig. 2, j(k) 30, is 
5 the complex valued time discrete base band representation of the received signal 16, 18. That 
is 



10 



25. 



0) y(A0 = rt(ft) + ir Q (A0 A = 1,2,... 



describes the possibly vector valued signal y. Here, i is the imaginary unit 
defined by i 2 = -1. The quantity y(Jt) 30 is often referred to as the multi channel time discrete, 

« 

digital base band communication signal. Examples of such modulation schemes are BPSK 

15 (Binary Phase Shift Keying) and one-dimensional amplitude modulation. If the transmitted 

« 

modulation form is MSK or OMSK, then the input signals to the Derotation block are 
derotated. The derotation works as follows. First, form thetjyo possibly vector valued 

— — y 

sequences ^ ^ 

/' 

20 . (2) y(fc) = ri(k) + ir Q (l:) *-l,2>... 

. - 

is the complex valued sequence constructed. Then, the derotated sequence 



(3) y(*0 =M~*y(fc), fc = i,2,... - z _i 

which is the output from the Derotation block 20 in the second case, is constructed. 

The output of the Derotation block 20 in Fig. 2 is die input to a Separation block 
22. The Separation block 22 separates the input signal in its real and imaginary parts: 



30 (4 ' 5) x r (*)=Re{y(*)}. 

x .(fc)=Im{y( fc )> 
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Here, Re{.} and Im{.} denote the real and imaginary parts of the indicated quantity, 
respectively. The vector, x(k) 32, formed from the real and imaginary parts of y(k) 30 is the 
output of the Separation block 22: 



(6) 



*«> - [ $ ] ■ 



The vector valued sequence x(k) 32 contains the transmitted signal, co-channel interfering 
users, additive noise and other disturbances. 

An Estimator block 2$ in Fig. 2, receives the outputs of the Separation 22, Filter 

* 

10 24, andDetector 26 blocks as inputs. Said Estimator block 28 performs one or more of the 
following: fine synchronisation, fine frequency ofiset estimation/compensation, data model 
estimation, estimation of channel and interference parameters. One or more of these quantities 
are outputted to the Derotation 20, Filter 24, and Detector 26 blocks. The estimation procedure 
uses knowledge about redundancy in the transmitted symbols to determine the unknown 

15 quantities. Redundancy in the transmitted symbol sequence can be in the form of a training 

» 

sequence, a preamble, synchronisation symbols, known payload, etc. There are several 
embodiments possible of the Estimator and Detector blocks as described below. A Filtering 
block 24 receives as input the real valued vector, x(k) 32, from the Separation block 22, and 
the coefficients in the data model (h, W, Q) from the Estimator block. 
20 The Filtering block 24 performs a multi channel whitening operation. This is not 

a traditional whitening filter as it takes into account to second order properties of the noise and 
co-channel interference. An output of the Filtering block 24 is input to the Detector 26 and 

• " * * * 

Estimator 28 blocks. 

A multi channel linear receiver that achieves interference rejection is described 

25 below. For sequence estimation, the time-discrete channel to the user of interest is modelled as 
an FIR-filter. Interference and noise are together modelled as a vector-valued autoregressive 
(VAR) system, see D. AsztSly and B. Ottersten, "MLSE and spatio-temporal interference 
rejection combining with antenna arrays" , in IX European Signal Processing Conference 
(EUSDPCP -98), Island of Rhodes, Greece, September 8-11 1998 and R. A. litis, "A GLRT- 

30 based spread-spectrum receiver for joint channel estimation and interference suppression" , 

r 

IEEE Transactions on Communications, 37(3): 277-288, March, 1989. The importance of 
telling the interference as temporally coloured, and non-circular, has been noted in, D. T. 



nun 
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M. Slock, "An interference cancelling multichannel matched filter", in Communications 
Theory Mini-Conference in conjunction with Globecom, pages 309^312, Paris, April 1997, 
IEEE. 

After whitening the received data samples, maximum-likelihood sequence 
detection, or reduced complexity sequence estimation, is performed by applying a Viterbi 
algorithm to the whitened sequence. 

* 

Novelty j rf the approach in accordancewith the present invention lie s in e fficie nt 
combiMdon of the recei ver structure and the exploitation of the signal constellation. A gain 
with the approach compared to maximum-likelihood detection without exploiting the signal 
constellation is significant. Methods and ideas presented herein are thus applicable to receivers 
that do not employ antenna arrays and where the modulation waveform is non-circular. A 



£ «.* T 1 ^, , 

y ~ l/'r 



GSM system, which is well known, is one such specific example. 

In contrast to methods exploiting antenna arrays the method according to the 
invention claimed needs no change of hardware. The present invention can of course also be 
used to further enhance the performance of receivers using antenna arrays and/or temporal 
oversampling with respect to the symbol period. To enhance the performance of several blocks 
in the radio receiver the modulation scheme in several systems for mobile communication, like 
GSM, can be exploited. In environments where interfering users are the dominant disturbance 
term largest performance improvement is obtained. 

In one embodiment, the multi channel Filter block 24 is a vector-valued linear 
IIR (Infinite Impulse Response) filter. In a preferred embodiment, the Filter block 24 is a 
multi-dimensional FIR filter with the Z-transform. 



JC 



(7) 



W(z) = I + £w(fc}z-*. 



The output of the filter is denoted with x'(*) 34. The filter performs the 
following operation: 



30 



K 

(8) x'(A) - x(fc) + Y, W(Qx(Jb - 1). 



Note that all the quantities above are real valued. 
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The Detector block receives as inputs x'(k) 34 from the Filter block 24 and the 
data model (h, W, Q) from the Estimator block 28 The outputs from the Detector block 26 
are estimates of the transmitted symbols. The Detector block 26 can have several 
embodiments. 

5 In one embodiment, taking a binary decision performs the detection; symbol by 

symbol based on a scalar input signal x*()t) 34, for example: 

x'(Jk) > a u r was transmitted 
(9, 10) -x!(k) < 0 a w 0" was transmitted 

In a preferred embodiment, a multi channel maximum likelihood sequence 
10 estimator (MLSE) is used for detection. It appropriately takes the data model, (h, W, Q) into 
account and produces as an estimate the most likely symbol. MLSE may be implemented with 
the well-known Viterbi algorithm or may be approximated with reduced complexity versions 
thereof. 

Below, embodiments for the Estimator 28, Filter 24, and Detector 26 blocks are 
15 described. To simplify the development, a data model for the received signal is first 

introduced. A common linear model for the communication channels is, in the discrete time 
domain, a linear FIR channel from the transmitter to the receiver. This modej for the vector 
and real valued signal x(Jt) is herein provided as 

M 

20 (11) *(*) - Yj h ® $ ( k - 0 + *(*). 

1=0 

where n(k) is the noise including receiver noise, model mismatch, and the 

interference from other users. Note that (11) is a real and vector valued description of a 

received communication signal. As is cleared out below, this description is more general and 

25 facilitates more interference suppression than-a^onoplex valued system with half the 

• . f \ \ . 

dimension. To detect the symbols of interest, s(k), the channel and the properties of the noise 

should preferably be known. Next, the noisei>rdjtfrties are described. 

The temporal correlation of noise term in (11) deserves special attention. Most 

often; this correlation is neglected. This is an appropriate model when no interference is 

30 present or when the delay-spread of the channel is small. However, when the delay spread is 

significant, the signal of the interfering users will, since they also satisfy an equation like (11), 

be temporally correlated. In one embodiment, the noise, n(Jk), is modelled as a vector-valued 
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moving average (MA)-process with extra noise input. Once the parameters of the MA-process 
are estimated, n(k) is whitened and the transmitted sequence is estimated using the Viterbi 
algorithm on the whitened received sequence. This detector has excellent performance, 
however, the drawback with such a procedure is that the number of states in the Viterbi 
5 equaliser will grow without bound as the number of transmitted symbols increases. For long 
block lengths this is not a practical approach. 

In the preferred embodiment, the interference and noise are modelled as a vector- 
valued autoregressive process of order J£, VAR(A), according to 

K 

io d2) *(*) + £ w - o = e (*)« 

where e(k) is a temporally white and Gaussian distributed process with spatial 

covaiiance 

15 (13) E{e(t)e T (t)} = Q. 

In one embodiment, the model order is determined using one o^ several well- 
known model order selection criteria. In an alternative embodiment, the model order is a 
receiver design parameter and is fixed to a prespecified value. In yet another embodiment, the 
20 model order is adapted during communication according to the bit error rate on the channel. In 
yet another embodiment, multi models are estimated and used in the filtering and detection 
stages, the model that gives the lowest bit error performance will produce the symbol 
estimates. 

Bit error performance can be evaluated by error detection and/or error correction 
25 codes, which are common in digital communication systems. 

* 

The above noise model is of course only an approximation of the real 
environment with noise and interference. The advantage with this particular choice of noise 
model is that the detection of the symbols is relatively straightforward. 

A maximum likelihood approach to estimate the transmitted symbols assumes 
30 knowledge of the parameters in the model. Here, it is outlined how to estimate the parameters 
in a system with a known symbol sequence, for example, a training sequence. During the 
training sequence both the channel and noise parameters are to be estimated. The assumptions 
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11 

made on the communication system is that the signal of interest and the interferes are burst 

■ • 

synchronised and that the scenario is time-invariant so that the parameters estimated during the 
training period are valid throughout the complete information burst. 

An example of such a system is a GSM system where the base stations are 
synchronised. It is of course possible to use the described receiver for the case of 

* 

unsynchronised base stations, some loss in performance will be experienced. It is also possible 
to extend the method to handle the case of unsynchronised base stations as well. However, in 
that case the VAR-model needs to be updated recursively during the information burst 
Adaptive estimation procedures for the VAR-model are known in the art. 

To outline the estimation procedure equations (26, 27) and (31) below are 
combined and the relation 



K 



(14) 



x(A) + £ W(Qx(Jb - 1) = Hs(Jfc) + n(fc). 



15 



20 



25 



30 



Note, that the relation is linear both in W(l),..., W(£) and H. Measured 
received sequence x(k) and the data symbols s(k) are both known quantities during the training 
sequence. This results in a simple least squares fit when simultaneously estimating the VAR- 
parameters in W(l),..., W(K) and the channel coefficients in BL 



(15, 16) 



*=1 



K 



x(A) + £W(Z)x(*-/)-Hs(*) 



where N is the number of samples for which the minimisation is performed and 
W denotes die model parameters W(l),..., W(/£). An estimate of Q is given by the sample 
covariance of the residual, ^(4). 

a?) Q=4x>(*)* r (*) 

Jb=l 

In most systems the transmitter and receiver frequencies differ slightly. Since this 
disturbs the receiver, this frequency offset must be estimated and compensated for. It is 
assumed that the signals, r,(Jk) and r Q (Jt) have a coarse frequency offset compensation according 
to what is known in the art Fine frequency offset estimation can be performed in the Estimator 
block. Output from the Filter block 24 will have reduced interference and t£us, this method 



a 
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will perform better than methods that do not take interference into account. Similarly to (6) the 
vector x'(fc). is partitioned and denoted: 



(18) : *w - [ iU 1 • 



Similarly, the following notation is introduced corresponding to the partitioned. 



equation (14): 



< 19 > Rw] - *<*> + E wow* - o = g$] + m 

* 

* • 

ft 

An estimate of the frequency offset, w 0 , may now be obtained by minimising 



(20) 



N' 

&o = argminV" 



[xf(fc)l _ [H r S(*) cos(o; 0 Ar)* 
W(*)J [H |S (A:)sin(cj 0 Jfe) 



where N' is the batch of data over which the minimisation is performed. In one 

* « 

lbodiment, the channel matrices H, and H, in (20) are partitioned from the estimated channel 
matrix obtained from the channel estimation procedure, the filtering matrices W are obtained 



Him 



from the channel estimation procedure, and the estimated symbols are the oujput of the 
Detector block. 

In an alternative embodiment, the minimisation is performed using the training 
data. In this case the minimisation is performed simultaneously over W, Hp Hj and o Q 



(21) min 

W,H ri H if u/ 0 



E 

*=1 



As an alternative, the frequency-of&et estimation is performed on training data 
assuming a trivial channel and filter (identity operation). In yet another embodiment, the 

■ * 

frequency-offset estimate is adaptively estimated using well-known adaptive versions of the 
criteria above. In general the adaptation rate of the frequency-offcet estimate is slower than that 

of the data model parameters. The frequency offset can be compensated for in the Derptation 

* 

block 20 by letting 

(22) y(k) = (rt(k) + ir q (k))e- iu ' k k =1,2,... 
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Synchronisation or timing estimation must be performed to be able to perform 
data detection. It is assumed that the signals, r^k) and r Q (Jfc) are coarsely synchronised, i.e., the 
location of the training data is roughly known. This coarse synchronisation is performed as is 
well known in die art. 

5 A fine timing estimation can be performed in die Estimation block. Output from 

the Filter block will have reduced interference and thus, this method will perform better than 
methods that do not take interference into account. An estimate of the timing, r 0 , may now be 
obtained by minimising 



(23) 



f 0 « argmin V) ||x'(A, r 0 ) - Hs(*)f 



Jb=l 



20 



25 



30 



miiiiiiiiiii 



where fT is die batch of data over which the minimisation is performed. In one 
lent, the channel matrix H and the filtering matrices W in (23) are obtained from the 
channel estimation procedure with the coarse timing. The estimated symbols in (23) are the 
15 output of the Detector block 26. 

In an alternative embodiment of the present invention, the timing estimation is 
performed using the training data. In this case the minimisation is performed, simultaneously 
over W, H and Tq. 



N 



(24) 



ain \\ 

,H,to > • * 



mm 
w 



K 



x(A, to) + £ W{l)x(k - /, To) - Hs(Jfc) 



As an alternative, the timing estimation is performed on training data assuming a 
trivial channel and filter (identity operation). In yet another embodiment, the timing is 
adaptively estimated using well-known adaptive versions of the criteria above. In general the 
adaptation rate of die timing estimate is slower than that of the data model parameters. In this 
and die following section it is assumed that the coefficients in the data model (h, W, Q) are 
known/estimated and the notation A will be omitted. 

As mentioned earlier, in one embodiment, the detector consists of a temporal 
whitening filter followed by a MLSE constructed for a temporally white noise process. 

* * 

Because the noise and interference are modelled together as a VAR-process, the whitening 
filter is a multi-dimensional FIR filter with the Z-transform 
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(25) 



W(z)«I + £W(*)*-* 



To formulate the whitened process mathematically, denote the output of the 
whitening filter with x'(A), that is, 

K 

( 26 > x'(Jfe) = x(Jfc) + 53 W(0x(* - /). 

The output of the whitening filter satisfies the equation 



10 (27) x'(fc) = £ - ') + 

1=0 

where e(Jt) is the process noise in the VAR-model of the disturbance. The new 
channel, IT, is of order M + Kznd is the convolution of the original channel with the VAR- 
model. Hence, the new channel IT is 

K 

(28) ti(k) = W(i)h(& - J), k « 0, . . . , M + if, 

i==0 



where W(0) = I. Note that the new channel is a finite order mdti-dimensiDnal 
HR-filter. This concludes the discussion of the whitening of the received data. 
20 To formulate the MLSE we construct the 2 x (M + K) channel matrix H 

according to 

a9) H'=[h'(0) ... b!{M + K)] 

25 and the symbol vector 

(30). s(*) =[*(*) ... s{k-M-K) ] r . 



With these notations the whitened two-dimensional system in (27) is compacdy 



30 written 



(31) *(Jfc) « H's(fc) + e(Jb) 



15 



Finally, the maximum likelihood estimate of the received sequence is 

(32) {s(k)} = arg min 

{*(*)} 

EWO " H's(/)) r Q-i(x'(0 - H's(/)). 

Summation is done for all the symbols s(i) to detect and also including parts of 
die training data (if available) to provide initial and/or final conditions. Minimisation is 
preferably performed with the Viterbi algorithm, which is well known in the art. 

Outputs from the detector, s(k) are either hard or soft estimates of die input data 
10 symbols^. These hard or soft decisions are used in die deinterleaving, error 
correction/detection before die final information symbols are estimated. 



t 
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Claims 

1 . A method for a digital receiver exploiting second order statistics for adaptive co- 
channel interference rejection in wireless communication, having digitally I, in phase, and Q, 
quadrature, branches of a received transmitted signal as input to the receiver, a coarse 

5 synchronization and a coarse frequency offset compensation have being performed on said 
signal, and comprising a means for denotation, means for separation, means for filtering, 
means for estimating and means for detecting transmitted symbols in the received signal, 
characterized in that it comprises the following steps: 

derotating the signal in said means for denotation, outputting a complex valued time 
1 0 discrete base band representation of the received signal if a used modulation form is one 
dimensional, else if a used modulation form is multiple dimensional then construct, from 
possibly two vector valued signal sequences, a complex valued signal sequence and 
outputting said complex valued signal; 

; separating said output signal in said means for separating it into its real and imaginary 
1 5 parts which vector is output from said means for separation containing the transmitted signal, co- 
channel interfering user signals, additive noise and other possible disturbance;, 

estimating in said means for estimation, using the outputs from said means for 
separation, means for filtering, and from said means for detecting as input, and performing on 
one or more of the following quantities: fine synchronisation, fine frequency oflset estimation or 
20 compensation, or data model estimation, outputting one or more of said quantities to said means 
. for derotating, filtering, and detecting; 

filtering in said means for filtering, using the real value output from said means for 
separation and coefficients in said data model output from the means for estimation and 
performing a whitening operation on the signal additionally taking into account a second older 
25 properties of noise and co-channel interference, the filtered signal being output to the means for 

■ * 

estimation and detection; and 

detecting transmitted symbols in said means for detecting, making use of the output 
from said means for filtering and the data model from the means for estimation, thereby 
improving co-channel rejection in wireless communication, thus making it possible to increase 
30 the number of channels for frequencies used. 

2. A method according to claim 1 characterized in that said filtering means 
perform a whitening operation that takes into account the second order statistical properties of 
the noise and co-channel interference. 
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3. A method according to claim 1-2 characterized in that detection is 
performed by taking a binary decision, symbol by symbol based on a scalar input signal. 

4. A method according to claim 1-2 characterized in that linear 
equalization is performed followed by symbol by symbol detection. 

5 5. A method according to claim 1-2 characterized in sequence estimation, 

is performed on the whitened sequence. 

6. A method according to claim 5characterized in that the sequence 
estimation is performed by using a filtered channel and a residual covariance matrix. 

7. A method according to claim 1-2 characterized in that detection is 

1 0 performed by a reduced complexity sequence estimation performed on (he whitened sequence. 

8. Amethodaccordingtociaiml-7characterized in that a model for a 
communication channel is a linear FIR channel. 

9. A method according to claim 1-8 characterized in that amodel for 
interference and noise is a linear vector valued moving average process. 

15 10. A method according to claim 1-8 characterized in that amodel for 

interference and noise is a vector-valued autoregressive process. 

11. A method according to claim 1-10 ch ar act erize d in that data model 
parameters are estimated using training data. 

12. A method according to claim l-10characterized in that data model 
20 parameters are estimated using estimated symbols. 

13. A method according to claim 1-12 characterized in that the model order 
is prespecified. 

14. A method according to claim 1-12 char acterized in that the model order 
is chosen according to a model order selection criteria. 

25 15. A method according to claim l-12ch aracter ized [ in that the model order 

is adapted during communication according to error rate on the channel. 

16. A m^hod according to claim l-12character ized in that multi models 
are estimated and used in the filtering and detection means, the model resulting in the lowest 
error rate providing the symbol estimates. 
30 17. A method according to claim 1-16 character ized in that timing 

estimation is performed on training data. 

18. A method according to claim 1-16 c h ar a ct e r ized in that timing 
estimation is performed using estimated symbols. 
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19. A method according to claim 1-16 characterized in that timing 
estimation is performed simultaneously as the data model parameters are estimated using training 
data. 

20. A method according to claim l-19characterized in that frequency offset 
5 estimation is performed on training data. 

21. A method according to claim 1-19 c h a r a c t e r i z e d in that frequency offset 
estimation is performed using estimated symbols. 

22. A method according to claim 1-19 characterized in that frequency offset 
estimation is perfoniied simultaneously as the data model parameters are estimated using training 

10 data. 

23. A receiver exploiting second order statistics for adaptive co-channel interference 
rejection in wireless communication, having digitally I, in phase, and Q, quadrature, branches 
of a received transmitted signal as input to the receiver, a coarse synchronisation and a coarse 
frequency offset compensation have being performed on said signal, and comprising a means 

15 for derotation, means for separation, means for filtering, means for estimating and means for 
detecting transmitted symbols in the received signal, characterized, inthatit 

* 

comprises: 

derotating means for derotating the signal, and outputting a complex ^alued time 
discrete base band representation of die received signal if a used modulation form is one 
20 dimensional, else if a used modulation form is multiple dimensional then construct, from 

possibly two vector valued signal sequences, a complex valued signal sequence and outputting 

■ 

said complex valued signal; 

separating means for separating the derotated output signal into its real and imaginary 
parts which vector is output from said means for separation containing the transmitted signal, co- 
25 channel interfering user signals, additive noise and other possible disturbance; 

estimating means using the outputs from said means for separation, means for filtering, 
and from said means for detecting as input, and performing on one or more of the following 
quantities: fine synchronisation, fine frequency offset estimation or compensation, or data model 
estimation, outputting one or more of said quantities to said means for derotating, filtering, and 
30 detecting; 

filtering means using the real value output from said means for separation and 
coefficients in said data model output from the means for estimation and performing a whitening 
operation on the signal additionally taking into account a second order properties of noise and 
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co-channel interference the filtered signal being output to the means for estimation and 
detection; and 

detecting means for transmitted symbols in the signal, making use of the output from 
said means for filtering and the data model from the means for estimation, thereby improving co- 
5 channel rejection in wireless communication, thus making it possible to increase the number of 
communication channels for frequencies used 

24. A receiver according to claim 23 characterized in that said filtering 
means perform a whitening operation that takes into account the second order statistical 
properties of the noise and co-channel interference. 
10 25. A receiver according to claim 23-24 characterized in that detection is 

performed by taking a binary decision, symbol by symbol based on a scalar input signal. 

26. A receiver according to claim 23-24 characterized in that linear 
equalization is performed followed by symbol by symbol detection. 

27. Areceiveraccordingtoclaim23-24character ized insequence 
1 5 estimation, is performed on the whitened sequence. 

28. A receive according to claim 27 characterized in that the sequence 
estimation is performed by using a filtered channel and a residual covariance matrix. 

29. A receiver according to claim 23-24 characterized in |bat detection is 
performed by a reduced complexity sequence estimation performed on the whitened sequence. 

20 30. Areceiveraccordingtoclaim23-29characterized in that a model for a 

communication channel is a linear FIR channel. 

m 

31. Areceiveraccoidingtoclaim23-30characterized in that a model for 
interference and noise is a linear vector valued moving average process. 

32. Areceiver according to claim 23-30 characterized in that a model for 
25 interference and noise is a vector-valued autoregressive process. 

33. A receiver according to claim 23-32 char acterized in that data model 
parameters are estimated using training data. 

34. A receiver according to claim 23-32 c h a r a c t e r i z e d in that data model 

• * 

parameters are estimated using estimated symbols. 

30 35. A receiver according to claim 23-34 characterized in that the model 
order is prespecified. 

36. A receiver according to claim 23-34 characterized in that the model 
order is chosen according to a model order selection criteria. 
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37. A receiver according to claim 23-34 characterized in that the model 
order is adapted during communication according to error rate on the channel. 

38. A receiver according to claim 23-34 characterized in that multi models 
are estimated and used in the filtering and detection means, the model resulting in the lowest 

5 error rate providing the symbol estimates. 

39. A receiver according to claim23-38 characterized in that timing 
. estimation is performed on training data. 

40. A receiver according to claim 23-3 8characterized in that timing 
estimation is performed using estimated symbols. 

10 41 . A receiver according to claim 23-3 8characterized in that timing 

estimation is performed simultaneously as the data model parameters are estimated using training 
data. 

42. A receiver according to claim 23-41characterized in that frequency 
offset estimation is performed on training data. 
15 43. A receiver according to claim 23-41 characterized in that frequency 

offset estimation is performed using estimated symbols. 

44. A receiver according to claim 23-41 characterized in that frequency 
offset estimation is performed simultaneously as the data model parameters arjb estimated using 
training data. 

20 
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